Scaling up the extrinsic curvature in gravitational initial data 
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Vacuum solutions to the Einstein equations can be viewed as the interplay between the geometry 
and the gravitational wave energy content. The constraints on initial data reflect this interaction. 
We assume we are looking at cosmological solutions to the Einstein equations so we assume that 
the 3-space is compact, without boundary. In this article we investigate, using both analytic and 
numerical techniques, what happens when the extrinsic curvature is increased while the background 
geometry is held fixed. This is equivalent to trying to magnify the local gravitational wave kinetic 
energy on an unchanged background. We find that the physical intrinsic curvature does not blow 
up. Rather the local volume of space expands to accommodate this attempt to increase the kinetic 
energy. 

PACS numbers: 04.20.Cv 



I. INTRODUCTION 

Initial data for the Einstein equations is usually con- 
structed by the conformal method. One is given 'free' 
data and rescales it to get the physical initial data. This 
is necessary because the Einstein initial data has con- 
straints. A comprehensive discussion of the constraints 
can be found in [Ij , especially in Chapter VII. Interesting 
physics tends to occur at the boundaries of the space of 
free data: one gets at the very least some insight into the 
limitations of the conformal method. 

The free data consists of a 'base' metric, a Riemannian 
3-metric, a divergence- free, trace-free symmetric tensor ( 
a TT tensor), and a scalar which is the trace of the ex- 
trinsic curvature. If the domain is a compact manifold 
without boundary we often choose the scalar as a con- 
stant. We then use a conformal transformation to solve 
the constraints. If neither the TT tensor nor the constant 
vanish, we can always find the appropriate conformal fac- 
tor ■ Parts of the boundary of the free data are easily 
acccsible. We can scale any one of the three parts by 
multiplying it by a constant and letting the constant be- 
come either large or small. What effect has this on the 
physics? For example, does the solution just cease to ex- 
ist, does the volume of the spacetimc blow up (or shrink 
to zero), do apparent horizons (which may be interpreted 
as cosmological horizons) appear? 

Maxwell's equations also have constrained initial data. 
By choosing the Maxwell free data as a pair of 3- vectors, 
(A,F), a parallel can be seen between the two fields. 
The magnetic and electric fields can be generated via 
B = V x A and E = F — VV, where V is a scalar 
function chosen to satisfy V 2 V — V;F\ It is clear that 
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if we multiply cither A or F by any constant, B or E will 
be multiplied by the same amount. As a result, this kind 
of rescaling can be used to increase the electromagnetic 
energy density (E 2 + B 2 ) without limit. 

Can we perform such a rescaling in General Relativity? 
Can we increase the gravitational wave energy at will by 
multiplying any part of the free data by a large constant? 
Because the gravitational wave energy interacts in a very 
nonlinear way with the geometry, it is not clear what 
happens with the Einstein equations. 

In this article, we discuss one such rescaling of the 
gravitational free data. The initial data consists of a 3- 
metric, gij, and a symmetric tensor, K lJ which is the 
extrinsic curvature of the slice. This is essentially the 
velocity of the 3-metric. K %3 is generated from the TT 
tensor and the constant. The constant is the Hubble 
constant, it represents a uniform expansion or contrac- 
tion of the space, while the TT tensor (which has only 
2 degrees of freedom per space point) can be interpreted 
as the gravitational wave velocity. We multiply the term 
that represents the gravitational wave velocity by a large 
number and track the solution. It turns out that the lo- 
cal energy density remains constant, while the volume of 
the universe expands. 

This article is a combination of analytic and numeri- 
cal works. Although we study a highly nonlinear elliptic 
equation, nevertheless we obtain some interesting ana- 
lytical results. In addition, we can accurately solve the 
equation numerically. This numerical work not only con- 
firms the analytic results, where available, but also indi- 
cates the behaviour of solutions in situations where we 
cannot prove anything. 

In Section II we confine ourselves to describing the con- 
formal method of solving the constraints. In Section III 
we show that if the extrinsic curvature vanishes nowhere 
on a compact manifold and we increase it, the conformal 
factor uniformly blows up. However, when the extrinsic 
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curvature vanishes somewhere, the analysis in Section III 
is no longer valid. To investigate this special case, we re- 
vert to a spherically symmetric toy model, deriving some 
analytic results in Section IV and showing the numeri- 
cal work in Section V. We distinguish between data sets 
where the extrinsic curvature vanishes in finite regions 
and when it only vanishes at an isolated point. Here we 
supply strong evidence that we do not get blow-up in re- 
gions of vanishing extrinsic curvature, and that we get 
slow blow-up when the extrinsic curvature vanishes at 
an isolated point. We conclude with a summary and an 
outline of future work. 



II. SOLVING THE EINSTEIN CONSTRAINTS 

Initial data for the Einstein equations consists of two 
parts: the first part is a manifold equipped with a Rie- 
mannian 3- metric gij, and the second is a symmetric 
tensor K ZJ on the same manifold. A' y is the extrinsic 
curvature of the 3-slice, i.e., the time derivative of the 
3-geometry. The metric and extrinsic curvature cannot 
be chosen arbitrarily. They must satisfy two constraints: 
the first is the Hamiltonian constraint, 



K lJ K,. 



K 2 



R 



K = 0, 



(1) 



where R is the 3-scalar curvature of and A is the 
trace of K l ° , i.e., A = g^K 13 . The second, the so-called 
momentum constraint, is 

ViK ij -WK = 0. (2) 

The terminology and notation comes from Q. 

The standard way of generating solutions to these 
equations is by means of a conformal transformation. On 
any given manifold, it is easy to construct TT tensors 
These are tensors that are both tracefree and diver- 
gencefree, gijK^ T = 0; ViA^? T = 0. Such TT tensors are 
conformally covariant. If one multiplies the given base 
metric by an arbitrary positive function (f> to con- 
struct a new metric , g^ = 4> 4 gij, then K^ T = <f>~ 10 Kj! T 
is TT with respect to g. Q 

Any K %3 which is the sum of a TT part and a constant 
trace, i.e., K lJ = K^ T + \Kg %3 , where A" is a constant, 
satisfies the momentum constraint Eq.Q. If we make a 
conformal tranformation of the metric, gij = (pgij, we 
find that K ij = <f)- 10 K^ T + \Kg^ again satisfies the 
momentum constraint with respect to the new metric. 
Note that K is not transformed; rather, it remains a 
given constant. In the conformal method, one uses the 
fact that we can freely choose the conformal factor to 
solve the Hamiltonian constraint. 

If we conformally transform the metric, g^ = <f> gij, 
we find that the scalar curvature transforms as 



R 



l R 



(3) 



We want the final metric and the final extrinsic curva- 
ture to satisfy the Hamiltonian constraint Eq.([T]), R — 



R-K^ T Kl T + |A 2 



0. This reduces 



to the famous Lichnerowicz-York equation [4|, |5[ 



v- 



R 



-K TT K ij < 



K 2 
12 V 
K 2 







(4) 
(5) 



where A 2 = Kj,' K'j', - 

This equation is very well behaved @]. In this article 
we focus on the situation where the topology of the 3- 
manifold is compact and without boundary. It can be 
shown that Eq. ([5]) has a unique positive solution if A ^ 
and if A^? T is not identically zero. It need not be non-zero 
everywhere; it is enough that it not vanish somewhere 0]. 
Let us remind the reader that A is a constant, while A^? T 
is a function. In the special cases, where either K = or 
K^ T = 0, we have an extra condition related to the sign 
of the scalar curvature. However, in the general case no 
such restriction applies. This existence result does not 
depend either on the metric, other than it be uniformly 
elliptic, or on the topology of the 3-manifold. 



To recapitulate 
free data 



i.9iji K-TTi 



(a K 13 



we start with a triplet, i.e., the 
A), and construct a new set 



K) 



9ij, 



~ 10 A"y T ,A) that satisfies the 



constraints. 



We would like to investigate the boundary of the set 
of free data. There exist parts of the boundary that are 
easy to reach. We can change the original triplet in a 
simple way by multiplying any one of them by a constant 
without touching the other two, because, for example, 
multiplying a TT tensor by a constant does not change its 
TT'ness with respect to a fixed metric. Of course, such a 
transformation will change the solution, </>, of Eq.([5]), and 
thus the data satisfying the constraints will be different. 

In this article we will consider such a rescaling. We 
pick one triplet (gij, A^? T , A), and use it to construct a 
family of free data of the form (gij,a 12 Kj! T ,K), where 
a is a running parameter. In particular, we wish to focus 
attention on the behaviour of <f> as a becomes large in 
order to see what happens to the physical initial data 
that emerges. This is a particularly interesting part of 
the boundary. 

One can regard A^ T as the velocity of the gravita- 
tional waves coded into the initial data. This rescaling 
seems to push up the gravitational wave energy content 
of the free data. Therefore, we are trying to increase the 
gravitational wave energy while controlling the rest of the 
geometry as much as we can. 

In the special case where A^? T vanishes nowhere, we 
can show that the conformal factor uniformly blows up. 
We demonstrate this in the next section. 
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III. A HARNACK TYPE INEQUALITY FOR 
THE CONFORMAL FACTOR 



We want to look at Eq.© 
R , 1 



8 



8 



-A 



12 



= 0, 



well behaved, has a regular positive solution. Let us look 
at what happens at the maximum of <f>, which we shall 
assume occurs at a point r = r max . The first two terms 



in Eq.© will be negative at 



3 ~ c 



so we get 



> 0. 



(10) 



where A 2 = Kf^ Kj! T - We are interested in the situation 

where we scale K^ T K^ T by a constant a 12 . Therefore 
we look at 



1 K 2 
8 12 Y 



0. 



(6) 



and we wish to show that q 
a becomes large. We write 
becomes 



scales linearly with a as 
t> = (j)/a and then Eq.© 



^0 + ^(^-7 



12 v 



) = 0. 



(7) 



We wish to solve the family of equations on a compact 
manifold without boundary. It turns out that the sign of 
the scalar curvature plays a minor role in the behaviour 
of the solutions. We can always set the scalar curvature 
to a constant value because of the Yamabe theorem Q, 
which tells us that any Riemannian metric on a compact 
manifold can be conformally transformed to a metric of 
constant scalar curvature (this really only makes sense on 
a compact manifold). The key quantity is the Yamabe 
number 



Y = inf 



Jl(V6) 2 + lR6 2 ]dv 

[fe e dv] 1 / 3 



(8) 



where the infimum is taken over all smooth functions, 
9. The sign of the Yamabe number fixes the sign of the 
constant scalar curvature. Since Eq.© is conformally 
covariant, and since conformal transformations form a 
group under composition, we can set R to a constant 
value without losing any generality. However, we do need 
to handle the three separate cases, Y > 0, R > 0, Y < 
0,R. < 0, and Y = 0, R = independently. In each 
case we will set the value of K 2 = 9. This choice does 
not change in any fundamental way the behaviour of the 
solution. 



This becomes 



[max^r < -A 



< - max A 2 . 
~ 6 



(11) 



The second inequality is needed because the location of 
fmax may well depend on a and, in general, does not co- 
incide with the maximum of A 2 . However, Eq.tfTTj) gives 
us a uniform upper bound on 4> as a constant independent 
of a. 

We get a uniform lower bound by looking at Eq.© 
when is a minimum, which we shall assume occurs at 
r = r min . At r = r mm we get 



- min A 2 <-A 2 \ 
6 ~ 6 1 



< 4^ 



+ [min0] i2 . (12) 



Consider the cubic equation 



x 1 

x 3 + 4— min A 2 = 0. 

or 6 



(13) 



This will have a positive root that increases as a in- 
creases. Fix a, a = ao, and find the positive root of 
Eq. (|13j) for a = ay. This number is a lower bound for 
[min0] 4 that is independent of a for all a > ao. In turn, 
this means that mine/) diverges at least as fast as a as 
a becomes large. Using the bounds on both min <f> and 
max0, we have shown that there exists a universal con- 
stant Cq independent of a such that 



max ( 



>C > 



min A 

max A 2 



1/3 



(14) 



The maximum and minimum of <j) both increase together 
proportional to a so that their ratio remains bounded 
independent of a. This can be regarded as a version 
of the Harnack inequality for the non-linear equation 
Eq.©. 



A. Y > 0, R > 

We assume that we are in the positive Yamabe class 
and set the scalar curvature R = +24; the specific num- 
ber can be chosen freely. Now Eq.((7]) reduces to 

V 2 0-30 + a 4 (-U 2 0- 7 -^ 5 )=O. (9) 
8 4 

Eq. © , because it is just a rescaled version of the original 
Lichnerowicz-York equation Eq.©, which is extremely 



B. Y<0, R<0 

We assume that we are in the negative Yamabe class, 
and set the scalar curvature R = —24. Now Eq.© re- 
duces to 

V 2 4> + 3j> + a 4 (lA 2 ql- 7 -h 5 ) = 0. (15) 
8 4 

Let us look again at what happens at the maximum of <fi. 
We again assume that this occurs at a point r = r max . 
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The first term in Eq. (jl5j) 

r max so we get 



V 2 0, will be negative at r 



and 



IaK 



if 



3(6 



This becomes 



~ 12 4[max( 



> o. 



< - max A 2 . 
6 



Consider the cubic equation 



4a; 2 1 , 9 
— ; max A = 0. 



or 



6 



(16) 



(17) 



(18) 



This equation has a positive root which is an upper bound 
for [max</>] 4 for any given a. This root decreases with in- 
creasing a. An easy way to confirm this is to differentiate 
Eq. (fT8"| with respect to a. We get 



8x . dx ^ 16ar 
a 4 da a 5 



0. 



(19) 



Using Eq.(T7]), we get 3a; 2 - 8x/a 4 = Ax /a 4 + 
max^4 2 /6a; > 0. This means that dx/da < 0, as re- 
quired. As in the positive Yamabe case, pick a value 
for a, a = ao, find the root of Eq. (fT5]) . and it will be a 
uniform upper bound of max^> for all a > ao. 

We get a uniform lower bound by looking at Eq. (fT5)) 
at the point where is a minimum. We assume that this 
occurs at r = r m i n . The first two terms in Eq. ([15[) are 
positive at the minimum. Therefore we get 



-minA 2 < -A 2 \ rmin < [min</>] 12 . 



(20) 



This inequality Eq. (|20[) gives the desired lower bound for 
[mine/)], which is independent of a. In turn, this means 
that min </> diverges at least as fast as a, as a becomes 
large. Using the bounds on both min0 and max0, we 
have shown that there again exists a universal constant 
Co independent of a such that 



-,1/3 



:A 2 



(21) 



and we again recover a Harnack inequality, but now for 
Eq.GSJ. 



C. Y = 0, R = 
Now Eq.([7]) can be reduced to 

W 2 4> + a 4 ( 1 -A 2 ^-^) = 0. 



(22) 



Showing the existence of uniform bounds in this case is 
even easier. We immediately get 



[max^] 12 < ^A 2 | rmax < imax^4 2 . 



(23) 



[minefil 12 > \a\ . >\m.inA 2 . (24) 
6 6 



Again we get 



:A 2 



1/3 



(25) 



D. Discussion 

The Harnack inequalities that have been derived in 
the last three subsections are clearly only valid when 
min(A 2 ) = mm(K^ T K^ T ) ^ 0. Since TT tensors are 
usually constructed by a decomposition method, one 
might think that it would be difficult to find TT tensors 
which vanish cither at points or in regions. Interestingly, 
we can construct such TT tensors [1(| and they cannot 
be ignored. It has not yet been possible to derive general 
results for such special TT tensors. Therefore we revert 
to considering only a spherically symmetric toy model. It 
turns out, in the spherical situation, that the case where 
A 2 vanishes in a region is easier to analyse than the case 
where A 2 vanishes at a point. When A 2 vanishes in a re- 
gion, we can prove that the minimum of cf> has an upper 
bound and does not scale with a. In the region where 
A 2 7^ we get the standard linear scaling with a. We use 
a combination of analytical and numerical techniques in 
the next two sections deal with this spherical model. 



IV. A SPHERICAL TOY MODEL: 
ANALYTICAL RESULTS 

We will restrict our attention to the case where the 
base 3-metric is spherically symmetric, and we replace 
the position-dependent function A 2 by an arbitrary 
spherically symmetric function because we know there 
exists no regular spherically symmetric TT tensor on 
flat space or on a spherical compact manifold without 
boundary 0. For convenience, we set K = 3. We start 
with a round 3-sphere of constant scalar curvature i?o (a 
natural choice is Rq = 24) and seek solutions to 



V 2 0-ii? o 



-a 12 A 2 



2K 2 
~2T 







1 



3</>+ -a 12 A 2 - 



-^ = 



(26) 
(27) 



for varying a. It is clear that when either a -> or 
a — > oo strange things happen. In keeping with the focus 
of this article, we will only consider here the case where 
a — > oo. 

Looking at Eq. (f2"Tj) , it is clear that (j> cannot remain 
regular as a — > oo, because the third term would diverge 
while the others remain regular. The results derived in 
Section III remain valid so we know that 4> is linearly 
proportional to a as long as A 2 is nowhere zero. We 
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need to deal with the situation where A 2 vanishes some- 
where. In particular, let us consider the situation where 
A 2 vanishes in a spherical region around the north pole. 

When dealing with spherical symmetry, we are free 
to take advantage of the fact that a round 3-sphere can 
be decompactified to flat 3-space, and that conformal 
transformations form a group under composition. In this 
picture the equation we wish to study is 



1 



a 12 A 2 c 



0, 



(28) 



with (/>—)■ at infinity. 

The conformal factor that maps flat space into a 
round sphere of scalar curvature equalling 24 is 9 = 



for any b > 0. The mapping that brings us 



from Eq. (|2"Tl) to Eq. (|2"5|) requires that A 2 and 4> be trans- 
formed. The transformation is that A 2 = B~~ 12 A 2 and 
4> = 6~ x 4>. In other words, <fi will be finite at the 'point 
at infinity' in the compact manifold ,while 4> ss Vb/r at 
the corresponding infinity in R 3 . 

A spherical region around the north pole corresponds 
to a disc < r < r\, where r is the standard confor- 
mally flat coordinate radius, on which A 2 = 0. It is now 
convenient to switch to the asymptotically flat picture. 
The equation we consider is Eq. (|28]l . and, in the disc 
< r < 7*i it reduces to 



value of <f> at the minimum. This uniquely determines the 
solution. The solution is 'U' shaped, blowing up twice at 
rj\ and rg. The bigger the min</>, the narrower the 'U', 
i.e., {tb — r^)(min^>) 2 is bounded However, we know 
that {ta < r±) and (rg > r 2 ) because the blow-up can- 
not occur in the range of validity of Eq (f2T)|) . Therefore 
tb — ta > r 2 — T\ , and so the value of min </> is bounded 
above. 

The more interesting case is when A 2 vanishes at a 
point rather than in a region. To repeat: we can show 
that when A 2 vanishes nowhere, we have uniform blow- 
up over the entire domain; when A 2 vanishes in a region, 
the minimum saturates. We conjecture that we will have 
behaviour which is 'halfway' between the two situations 
dealt with above. More precisely, we conjecture that the 
minimum will blow up with a, but at a rate which is 
slower than linear. The details will depend on the rate 
that A 2 goes to zero at the point. We have done some 
numerical modeling, both to confirm the analytic results 
and to investigate those situations where we can prove 
nothing concrete. These models will be discussed in the 
next section. 



A SPHERICAL TOY MODEL: NUMERICAL 
RESULTS 



(29) 



We can write down the solution of this equation explic- 
itly. It is 



(30) 



where a is a parameter. These are the functions that map 
one from flat space to the round hypcrboloid of constant 
negative scalar curvature. These functions blow up at 
r = a, and since we know that the total solution is reg- 
ular, we know that the blow-up must occur outside the 
range of validity of these functions. While we do not 
know a priori the value of a, and it will change with a, 
we do know that we have a lower bound for a, i.e., a > r\. 
The minimum value of <p for these special solutions occurs 
at the origin. There we get 



-4> = 4>{r = 0) 



< 



(31) 



Therefore min (f> does not blow up like a but reaches some 
limit, while max</> becomes unboundedly large. A similar 
argument holds when A 2 vanishes near the south pole. 

We can repeat this argument when A 2 vanishes on 
some belt t\ < r < r 2 . The equation, Eq. (f2T)|) . is still 
the same, but now is valid for t\ < r < r%. We cannot 
write down a set of explicit solutions, but we have much 
information about the solutions. Fix the location of the 
minimum, in this case in the interval (n, r 2 ), and fix the 



In this section, the equation we deal with is a one- 
dimensional elliptic equation with a highly nonlinear 
source term, 



V 2 <t) + S{4>) = 0. 



(32) 



We use a one-dimensional pseudo spectral method to 
solve this equation in two computational domains [l2| - 
[T3 |. We have a coordinate y <G [0, 2) replacing r € [0, oo). 
In the interior domain we have y = r,y G [0,l],r 6 [0, 1], 
while in the exterior domain we have y = 2 — 1 jr. y G 
(l,2),r G (l,oo). Therefore, infinity can be included 
when we put the computational grid at the point y = 2. 

We use Chebyshev polynomials as basis functions and 
the collocation points are 



r i 

Vj = l + cos[^-] in [1,2] 



(33) 
(34) 



where i = 0,1. ..n and j = 0, l,2...n. The two domains 
meet at i = on the interior domain and j = n on 
the exterior domain. The grid point j = touches the 
infinity of coordinate r. 

We require the solution to be C 1 . This means that 
4> and its normal derivative must match at the interface 
between the domains. The discrete equations are solved 
by the Newton-Raphson method. The resolution is taken 
as 100 on each domain. 
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A. A = 1 

We here consider the simplest case, where we choose 
the base metric to be a round sphere with scalar curva- 
ture equalling 24, and pick A 2 to be a global constant 
equalling 1. Then Eq. ([27|) 



1 



- 30 + -a 12 A 2 < 

8 

reduces to an algebraic equation 

, a 1 



<f = 0, 



i 12 



or in terms of the normalised 



a 12 =0, 



la, 



U2 



+ 



l S 1 

--6=°' 



(35) 



(36) 



— > i from below, or <f> — > 



As a — > oo, we expect </> 12 
0.8612992. 

The round base metric can be written as 



1 



9ij 



(1 + r 2 ) 2 13 ' 



(37) 



The computation is all done in the flat space, so we have 
to solve 



1 



a 12 A 2 c 



0, 



(38) 



where A 2 = 



(l+r a ) e ' 



At the risk of confusion we intro- 



duce a normalized <p = <p/a, then Eq. (|38|l becomes 



7 2 ±! 



1 



+ i-a 4 A 2 ^ 7 



(39) 



The relationship between <j> and ^ is <fi = vT + r 2 ^. 
We solve Eq. (f3"9")l numerically with A 2 



(l+r 2 ) 6 • 

We see that the solution is, as expected, of the form 
cf)' = C/\/l + r 2 with C depending on a. In Figure 1 
we present 4> on the compact manifold (by multiplying 
</>' by vT+r 2 ). We clearly see that we get a family of 
constant functions that asymptote to a fixed function as 
a becomes large. 

We list the values of C w.r.t. a as the table I. 



B. A vanishing at a single point 

A simple choice of A 2 that corresponds to A 2 vanishing 
at the south pole is to pick A 2 = r^rp with (3 > 6. Wc 
restrict our attentions to /3 = 7 and j3 = 10. We present 
(f> for several values of a with (3 = 7 and (3 = 10 in Figure 
2 and 3, where a — 1,2,3,4,5,10 gives rise to the red, 
green, blue, pink, light blue, and yellow lines respectively. 

It is clear that, in both cases, the maximum of <p, at 
the north pole, settles to a constant value independent 



0.9 

0.85 

0.8 
0.75 - 

0.7 - 
0.65 

0.6 
0.55 

0.5 



'1.000' 
'2.154' 
'4.642' 
'10.00' 
'21.54' 







Figure 1. ^ on S' 



a 


C 


1 


0.6680900040565467 


2.154434690 


0.8387341686204405 


4.641588834 


0.8601789420392659 


10 


0.8612470753749010 


21.54434690 


0.8612968154854488 


46.41972069 


0.8612991245657523 



Table I. It is a list of the scaling parameter a and the con- 
stancy of the standard solution (f> . It is clear that, as ex- 
pected, C asymptotes to 0.8612992. 



of a. This shows that the maximum of </> grows linearly 
with a. On the other hand, the minimum of </), at the 
south pole, where A 2 = 0, decreases with increasing a. 
Hence the minimum of <fi does not grow linearly with a. 

To analyse the behaviour of the minimum of <p as a 
function of a, we consider each of the two cases, i.e., 
(3 = 7 and (3 = 10. We conjecture that <p(9 = ir) scales 
with some power of a, for large a. We plot In 4>/ In a 
versus a for each of the two choices of f3. These are 
Figures 4 and 5. Each of the two curves flattens out 
for large a. These show that min (j) grows like a n with 



i.i 
1 

0.9 

o.s 
0.7 
0.6 
0.5 
0.4 



'1' 
'2' 
'3' 

T 
'5' 
10' 



Figure 2. 4> on 5* 3 with /3 = 7 
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Figure 3. cj> on S 3 with f3 = 10 



Figure 6. A 2 on S 3 



n ~ 0.71 when P = 7, and min</> grows like a n with 
n fn 0.405 when p = 10. 



growth rate seems to diminish smoothly as we move from 
P = 6 to p = oo. 
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Figure 4. In <f>/ In a vs a for ft = 7 




Figure 5. ln^i/lna vs a for j3 — 10 

If /3 = 6, A 2 does not vanish at the south pole and 
we get linear growth of <fi there. The limit p = oo corre- 
sponds to A 2 vanishing in a region near the south pole, 
and we get no growth at all. It is nice to sec that the 



C. A 2 vanishing in a patch 

We start by finding, on the round sphere, a smooth 
positive spherical function that vanishes in a region. We 
will use this as A 2 . We first find a cut function rj(x) 
which is defined on [0, 2] by 





209 - 2240a; + 10080a; 2 
-24640a; 3 + 35280a; 4 - 29568a; 5 
+ 13440a; 6 - 2560a; 7 
1 



x e [0,0.5) 



x e [0.5,1.0) 
x G [1.0, oo) 



This function smoothly interpolates between at a' = 
0.5 and 1 at x = 1 and is C 3 at each end. In terms of the 
rescaled coordinate function y which we introduced at 
the beginning of this section, we construct the following 
function 



f(y) 



7/(1 -y) ye [0,0.5) 

ye [0.5, 1.25) 
77(22/ -2) yG [1.25, 1.5) 

1 ye [1.5,2) 



This is a function that is 1 at the origin, goes to zero 
at y = r = 0.5, is zero in y G (0.5, 1.25), rises to 1 in 
y € (1-25, 1.5), and stays equal to unity out to y = 2, r = 
00. We choose A 2 = /(y)/(l + r 2 ) 6 . This corresponds to 
a smooth A 2 on the sphere, which equals 1 at both the 
north and south poles but equals zero in a central region. 
Figure 6 shows A 2 on the sphere. 

We solve the equation for a range of parameters a. 
We display the various solutions in Figure 7, showing the 
normalized <f> = <fi/a on S 3 . The parameters used are 
a = 1.000, 2.154, 4.642, 10.00, 21.54, 53.13, 79.37, and 
100.0. These correspond to the red, green, blue, pink, 
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light blue, yellow, dark blue, and mauve lines respectively 
on Figure 7. 



0.9 








Figure 7. <j> on S 3 

One can clearly see that the value of (f>, on the support 
of A 2 , tends to a stationary limit, while <f> collapses off 
the support of A 2 . This shows that <f> scales linearly with 
a, on the support of A 2 , while in the region where A 2 
is zero, <f> continues to diminish so that cj> approaches a 
stationary value. 

We wish to show that min (f> = a min (f> increases with a 
but approaches some fixed upper bound. The minimum 
occurs around 6 = 1.35. In Figure 8 we have plotted 
ln<f)/lna versus a at 8 = 1.35. One should compare this 
graph with Figures 4 and 5. This looks like a graph which 
is going to an asymptotic value of 0, which indicates that 
<fi, at 9 = 1.35, is heading for a fixed number, independent 
of a. 
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Figure 8. ln^/lna vs a 



VI. CONCLUSIONS 

We have shown that if we scale K^ T K TT ij by a 12 , we 
find that the conformal factor, in general, scales like a. 



However, the physical K^ T K^ T = c^- 12 a 12 K^ T K^ T = 
(f)^ 12 K^ T Kfj T , and as a becomes larger and larger <j> 

remains finite. This means that K^ T K^ T remains fi- 
nite. Hence the velocity part of the gravitational wave 
energy density remains bounded even though the cor- 
responding 'free' data blows up. On the other hand, 
the volume of space becomes unboundedly large, because 
Vs = <f> e '"Jd = o^i^y/g blows up. Therefore the total 
gravitational wave energy in a coordinate sphere becomes 
larger and larger while the local energy density remains 
bounded. 

We see that the gravitational wave energy inside any 
coordinate sphere increases like a e while the surface area 
increases like a 4 . This indicates that this family of initial 
data will eventually contain horizons. 

The numerics, when combined with the analytic cal- 
culations, show a coherent picture. If K^ T Kfj T has no 
zeros, then the conformal factor blows up uniformly. If 
the K^ T Kj 3 T is zero on a patch, we expect no blow up on 
this patch, but, nevertheless, we continue to get the stan- 
dard blowup elsewhere. If K^ T K^ T has an isolated zero, 
we will get blowup at this point, but at a rate slower than 
in the rest of the space. The rate of blowup is not uni- 
versal in this case, but depends on how quickly Kj, T Kf^ 
moves away from zero at that point. There seems to be 
a smooth transition between the 'slow blow-up' case and 
the 'no blow-up' case. 

There are a number of obvious extensions to this work. 
To find a real TT tensor, we need to abandon spherical 
symmetry and, at the very least, work with axially sym- 
metric data. We would expect that if we have a real 
K^ T K^ T , which either vanishes at a point or vanishes in 
a region, to get behaviour similar to the spherical model. 
We would be surprised if anything could be proven an- 
alytically; we are going to have to depend on numerical 
modeling. Work is in progress in this direction. 

It would be interesting to repeat this analysis in the 
asymptotically flat case. We would probably want to 
work with maximal initial data, i.e., K = 0, and just 
have a metric and a TT tensor as free data. It is clear 
that one can change the metric so that the ADM mass 
becomes unboundedly large and trapped surfaces appear 
[ll1 | . What happens if we blow-up the extrinsic curvature 
on a fixed background metric? Will we get the same 
behaviour? Preliminary investigations indicate that we 
do: the ADM mass diverges and trapped surfaces appear. 
We intend to investigate this further. 
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